Non local lagrangians: the pion. 
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We define a family of non local and chirally symmetric low energy lagrangians motivated by 
theoretical studies on Quantum Chromodynamics. These models lead to quark propagators with 
non trivial momentum dependencies. We define the formalism for two body bound states and 
apply it to the pion. We study the coupling of the photon and W bosons with special attention 
to the implementation of local gauge invariance. We calculate the pion decay constant recovering 
the Goldberger-Treiman and the Gell-Mann-Oakes-Renner relations. We recover a form of the axial 
current consistent with PCAC. Finally we study the pion form factor and we construct the operators 
involved in its parton distribution. 
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The strong interaction among hadrons is supposed to be described by Quantum Chromodynamics (QCD) 1] which 
is a field theory defined in terms of quark and gluon fields. While the asymptotic behavior of QCD is well understood 
and its proponents worthy of the highest recognition 0, the low energy behavior is still a subject of much scientific 
endeavor. Low energy physics seems to be ultimately governed by flavor dynamics. Confinement the property of 
QCD which describes how the dynamics based on color in the lagrangian transforms into a dynamics based on flavor 
for the physical states, and why these cannot exist with color charge, is still a subject of research and debate. This 
complex low energy behavior is described conventionally in terms of approximations to the theory, i.e., lattice QCD 
0, non relativistic QCD 0, 1/N C -expansion 6] or effective theories, i.e. Chiral Perturbation Theory 0, Heavy 
Quark Effective Theory Q, etc. Models turn out to be extremely useful in some instances when the other approaches 
are too complex, i.e, non-relativistic quark models bag models [Hill, Nambu-Jona Lasinio (NJL) model [3, 

C"| ■ chiral-soliton models [lil flH fl6l| . etc. Another method to study non perturbative physics in a lagrangian theory like 
Oh! QCD is to solve the Dyson-Schwinger equations. The application of this formalism to QCD becomes an enormously 
difficult task but progress, in understanding the theory from this appro ach, has been achieved The global 

<j_) \ color model [2(J , the extended non-local NJL model 0, |2^, |23, I24l25j , and models using separable interactions [2|J 
r"| , has been introduced as model realizations of QCD in a field theory formalism. 

One major problem is to understand the pion, because it is the system which contains all the ingredient of QCD: 
• i-h , asymptotic freedom, confinement and spontaneously broken chiral symmetry. Chiral symmetry governs the static 
properties of the theory, like the quark condensate, the mass and decay constant of the pion. The dynamics fixes the 
internal structure of the pion, which is accessible through the pion electromagnetic form factor. 

The pion form factor has been a subject of many studies. In relativistic quantum mechanics the pion form factor 



has a long history of debate |27], l28l l29l l3Ct l31| . One of the main problems is the choice of one of the various Dirac 
forms |32l 1331] . One of the most interesting conclusions is that results of a calculation depend not only on the Dirac 
form chosen but also on the frame chosen |3 ll l35| . The reason for this is that the truncation needed in relativistic 
quantum mechanics to calculate form factors breaks Poincare invariance. 

A way to avoid this problem is to work in a field theory formalism. The pion form factor has also been studied 
within Dyson-Schwinger Equations schemes by several authors |36l l37l l38| . The starting point is in most treatments 
the pion Bether-Salpeter amplitude calculated in the rainbow approximation. The pion form factor is calculated using 
the so called impulse approximation which considers only the triangle diagram, and the use of a dressed vertex for 
the photon. For the latter the Ball-Chiu j3(| expression for the vertex , or modified versions of it, have been used. 

The pion form factor has been studied also for the non-local NJL model . In this case the coupling of the photon 
is obtained by restoration of the gauge symmetry pH I24I lloj . 

Our aim is to construct a model for hadron structure and hadronic interactions developing a formalism which 
preserves the fundamental symmetries of the theory (chiral, Poincare and local electromagnetic gauge invariances) 
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and which incorporates information coming from fundamental studies of Quantum Chromodynamics. For this purpose 
we want to define a formalism which contains the physical intuition of model calculations and the lagrangian formalism 
of the effective theories. To achieve this we have found that the best suited scheme is to describe the physics by mean 
of a phcnomenological chirally invariant non local lagrangian. 

Working in a lagrangian theory, the two main ingredients in a non perturbative analysis involving the pion are: i) 
the quark propagator, obeying the Dyson equation; ii) the description of the pion as a bound state of a Bethe-Salpeter 
equation (BSE). Due to chiral symmetry the kernels of these two equations are not independent |4l| . Solving the 
Dyson equation for our lagrangian leads to momentum dependencies in the quark propagators through its mass and its 
wave function renormalization. In our scheme the gluons have been integrated out and we have only flavor interaction 
between quarks. Confinement is imposed by the structure of the quark propagator and by limiting the Fock space to 
color singlet states. The pion is obtained in a consistent way solving the BSE, and the Goldstone character of the 
pion is recovered. 

Our model can be seen as an extension of the non-local NJL model 0, |2^, 0, E3 > but with a particular philosophy. 
We consider the description of the quark propagator as the main ingredient. This is because the quark propagator is 
the first information that can be obtained from fundamental studies, as lattice QCD. Our lagrangian is the minimal 
extension which allows to incorporate the full momentum dependence of the quark propagator, through its mass and 
wave function renormalization. From this lagrangian we can explore what are the implications for other observables 
originated by changes in the quark propagator. 

Our formalism implements the coupling of the photon in a gauge invariant manner pH l24l l4p| . This allows to study 
the electromagnetic properties of the pion, which depend strongly on the quark-photon vertex. Usually this vertex is 
calculated by using the Ward-Takahashi identity 39] . This method fixes the longitudinal part of the vertex leaving 
the transverse part unconstrained. We show that this procedure does not guarantee local gauge invariance, while 
ours does. As an application we study the pion form factor, showing that the conventional impulse approximation, 
in which the form factor is calculated using the triangle diagram, is not consistent with local gauge symmetry. 

In models based on field theory formalism, the construction of the axial current and the definition of the pion decay 
constant need particular attention. In ref |43 a first expression for the axial current is given. In ref plll2^| additional 
contributions to the pion decay constant are included. In this paper we implement the coupling of quarks to the 
bosons in a gauge invariant manner following a procedure similar to the one used for photons. Then, we analyze the 
axial current and the pion decay constant. 

Our formalism is very effective for building operators describing observables in a consistent way. As an application 
we have studied the operators involved in the parton distribution of the pion. 

This paper is organized as follows. In section [H] we define our lagrangian, we discuss the quark propagator, and 
we fix the parameters in order to describe the adequate quark propagator obtained by more fundamental studies 
based on QCD. In section UTTl we describe the pion state. In section Hvl we study the quark-photon vertex and recover 
the Ball-Chiu ansatz but with additional contributions. In section we study the axial current and the pion decay 
constant. In section IvTI we apply the model to the study of the pion form factor. We show that a four quark-one 
photon vertex appears in a natural way. In section IVIII we obtain the contribution of this new term to the parton 
distribution operator. The last section contains the conclusions of our investigation. 



II. A PHENOMENOLOGICAL NON LOCAL LAGRANGIAN FOR HADRON STRUCTURE. 



Let us build a model which produces a non trivial momentum dependence in the quark propagator and preserves 
all the required symmetries: Poincare and chiral symmetry. This momentum dependence will arise from a lagrangian 
description and manifests itself as a quark momentum dependent mass and a quark momentum dependent wave 
function renormalization. Let us define the non local currents as 

Jo (x) = J d A yG (y) $\x+ \y^j Oi> (x - \y\ , (II.l) 
where the operator O is such that 

7°0 t 7° = O . {11.2) 

With these definitions the hermiticity of the currents ^ J a (x) = Jq (x)j implies that (— y) — G (y) . A local current 
corresponds to G (y) = S 4 (y) and thus the natural normalization for the functions G (y) is 

d i yG(y) = l. (11.3) 
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FIG. 1: Diagrammatic representation of the quark propagator. The non local four quark vertex is represented by a circle with 
structure. 

We build a lagrangian in terms of non local currents, preserving Ul (2) ® Ur (2) chiral symmetry, and producing 
the desired momentum dependences as 

C(x) = ij (x) (ip- mo) i> (x) + go [4 0) J s (x) + J5 ( x ) -h (x)] + (») (a?) , (H.4) 
where the currents are defined by 

Js (x) = J d 4 y Go (y) $ (x + ^y) i(,(x- \y^j , (II.5) 

■h (x) = J d 4 yG (y)i> (x + ^y^j ir^ip (x - , (II.6) 

J p (x) = J dS G p (y) $L+ i») (x - lyj , (II.7) 

where u d v — u (dv) — (du) v. The transformation properties of the non local currents are the same as those of the 
local ones. The first and second currents require the same Go (y) to guarantee chiral invariance. The third current 
is self-invariant under chiral transformations. The scalar current, Js, generates a momentum dependent mass, and 
the last current, the "momentum" current, J p , is responsible for the momentum dependence of the wave function 
renormalization. The pseudo-scalar current, J5, generates the pion pole. From now on, just for simplicity, we assume 
that all the G (y) functions are real. 1 

The interaction vertex obtained from the Lagrangian (|II.4|I automatically includes vertex form factors. Let us define 



G(p) = J d'ye^Giy) , (II.8) 

with the normalization condition 

G(p = 0) = l. (II.9) 
The full quark propagator is obtained from de Dyson equation, represented in Fig. ^ 

S ^ = 2 zWtx^ ' ( IL1 °) 

p - m + 2J (p) + le 

with 

X(p) = -a G ( P )-ta p G p (p) , (11.11) 

where the first term arises from the scalar current and the second from the momentum current. The constants ao 
and a p are directly related to the couplings go and g p , 

ao = 2go j ^ Go (p) Tr (iS (p))=i8 N c N f9o \ j-^ Go (p) ^ 2 - 7; , (H.12) 

(11.13) 



(2tt) 4 ^ " " " Jy J (2tt) 4 u ^V-m 2 (p)+i6 ' 

a p = 2 g p [ -^G p (p) Tr (iS (p) /}) = i 8 iV c iV /ffp / -^G p (p) 2 ^f^^. , 
J (2tt) 7 (2tt) P - m z (p) + le 



1 In a previous paper another derivative coupling, ( if) (x) ^id^ip (x) ) , was introduced l25l . This term is built with vector currents and 

J 2 
p- 



therefore, it produces different effects than our term J 2 In particular, it was used to reproduce vector meson dominance. 
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where Tr represents the trace in Dirac, color and flavor indices. For simplicity we work from now in the large N c 
limit. This is equivalent to the Hartree approximation which implies that only direct terms are taken into account. 

We can rewrite the momentum dependence of the quark propagator in a more standard way through a momentum 
dependence in the quark mass and in the quark wave function renormalization, 



$ + to (p) 

with 



s(p) = z(p) ^^r- • (H.14) 



m Q + a Q G Q (p) 

m(p) = - , , (11.15) 

1 - a p Gp (p) 

z (p)=i ■ ( IL16 ) 

1 - a p Gp (p) 

These relations between (Go (p) and G p (p)) and (to (p) and Z (p)) assures the self consistency of the solution of the 
Dyson equation. 

Eqs. (Ill - 1 51) and (|II.16|) show that the scalar current can give a mass to the quark even if the lagrangian contains 
no mass term, too = 0. This phenomenon is the spontaneous symmetry breaking mechanism which is similar to 
that taking place in the Nambu-Jona Lasinio model |13| . On the other hand, the momentum current gives rise to a 
momentum dependent wave function normalization. However, although it contributes to the mass, it is not able by 
itself to break spontaneously chiral symmetry. 

We shall be guided by fundamental studies of QCD and lattice parametrizations for building models for Go (p) 
and G p (p) . The natural way to proceed is to use the information coming from these studies to write ansatze for 
to (p) and Z (p). Then, transposing Eas. (|II.15|) and (|II.16(1 we obtain Go (p) and G p (p) . The values for ao and a p are 
determined from the normalization condition equation (|II.9|) 



= ~zjffj ~ m ° ' ( IL17 ) 

c^l-^y, (H.18) 

and, from Eqs. (|II.12(I and 13|) . we obtain the values for go and g p . 

These studies are performed in Euclidean space and therefore we will perform our calculations in this space. We 
use pe to represent the momentum in Euclidean space. 

Here Go (pe) , G p (pe) , Z (pe) and m (pe) are functions of p\. We impose that for p 2 E — ► oo, the mass goes to the 
current mass and the wave function renormalization to 1, 

m{p E ) — > m , (11.19) 

p 2 E ~*oo 

Z( PE ) 2 — 1. (11.20) 



Assuming that the integrals in Eqs. (jII.12|) and ljII.13|) are convergent, and looking at the behavior of the integrands 
for large values of p E we obtain that 

G {pe) — ► P E ° with a>2 + e (11.21) 
G p {p E ) — ► P E a with a>A + e. (11.22) 

Let us define Gq(pe) and G p (j>e) or alternatively m(pE) and Z (pe) ■ Much research has been carried out in 
the study of their functional shapes. We extract from these studies two well known scenarios based on different 
philosophies but equally consistent. 

The first scenario, which we will call SI, is based on the work of Dyakonov and Petrov |43(. They provide us with 
the momentum dependence of the quark mass term coming from an instanton model. They assume Z (pe) = 1 and 
work in the chiral limit (too = 0). Their results are well described by the expression 

A 2 X 3/2 



Pe 



m ( PE ) = too + a m ( A2 m ^ 2 ) , (11.23) 
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Case 


(qq) 1 '* (MeV) 


m (MeV) 


m-n (MeV) 


/x(MeV) 


< r 2 > (fm 2 ) 


SI 


-303. 


2.3 


137. 


82. 


0.41 (0.41) 


S2 


-285. 


3.0 


139. 


81. 


0.36 (0.40) 


Exp. 


-250. 300. 


1.5 ~ 8. 


135. ~ 140. 


92. 


0.44 



TABLE I: Results for < qq > 1 / 3 , mo, the corresponding pion mass, the pion decay constant and the mean square radius, 
< r 2 >, for the full vertices given by Eq. HI V.lOt and, between brakets, the Ball-Chiu ansatz for the two scenarios described 
in the main text. 



with A m = 0.767 GeV and a m = 0.343 GeV. 

The second scenario, which we call S2, corresponds to an alternative mass function obtained from lattice calculations 
as proposed by Bowman et al. 0, , 



to (p E ) = m + a r 



A 3 



(11.24) 



with A m = 0.719 GeV and a m — 0.302 GeV. In their lattice analysis the authors also look for the wave function 
renormalization constant. Their values are reasonably reproduced by 



Z (p E ) = 1 + a z 



A 



2 \ 5/2 



■Pe 



(11.25) 



with a z = -0.5 and A z = 1.183 GeV. 

In table [I] we show the values of some observables for the different scenarios. Among them, the quark condensate 
is defined by 



(qq) = ~i4N c 



d A p ( Z (p) to (p) 



to 



(2-7t) 4 \P 2 — m 2 (p) + ie p 2 — TOg + ie 
We stress that these values are obtained without any free parameter and therefore they are model predictions 



(11.26) 



III. THE PION MASS. 

In our formalism the Bcthc-Salpeter amplitude in the two body pion channel is defined as 

X 1 ( P , P)=iS{p+^p] i 75 r* fa (p) iS\p-^P 
where fa (p) is given by 

fa (p) - -i 2 5o Go (p) J -^4 G o ip') Tr (i 1 5 tH S (p' + ^ i j 5 r l fa (p') is(p'- \p 



(m.i) 



(III.2) 



which is represented in Fig. [5J Note that in equation ljIII.2f) there is no summation with respect to the isospin index 



The solution of equation I|III.2J| is straightforward and gives 

ip) = QnqqGo (p) . (111.3) 

The pion mass is obtained from the BSE, which can be easily rewritten in terms of the pseudo-scalar polarizability 

iV 
(2nf 



S^Ups (P 2 ) =-ij 0^G 2 Q (p') Tr (i 15 tH S (p' + l -P^j i l5 rU S (p' - Ip^j ) (III.4) 



and equation (|III.2I) becomes 



l = 2.g np S (P 2 



(III.5) 



FIG. 2: Diagramatic representation of the Bethe Salpeter equation. 



The normalization constant g^ qq is obtained by the usual normalization condition of the BSE which can be rewritten 

as 



1 



fdU 



PS 



Q 2 



\ dP 2 



(IIL6) 



P 2 —m 2 



As shown in tabled we obtain the physical pion mass for reasonable values of the current quark mass mo. 

The model realizes the Goldstone theorem. To see it explicitly we can go to the exact chiral limit, by choosing the 

reference frame where P^ = LM, Oj and taking the limit M — * 0. The explicit realization of the Goldstone theorem 

arises because chiral symmetry implies that we must use the same kernel in the Dyson equation for the mass, equation 
(|II.12ll . and in the BSE for the pion, equation l|III.2|) 41] . Notice that G p (p) is not constrained in the procedure. 



IV. THE QUARK-PHOTON VERTEX. 



In order to study the electromagnetic properties of the pion we describe the coupling of the dressed quarks to the 
photon. The usual approach to the quark photon vertex is to exploit the Ward-Takahashi Identity (WTI), which is 
a consequence of gauge invariance in QED. The WTI is satisfied order by order in perturbation theory and must be 
satisfied also non perturbativcly in order to have the right normalization for the quark-photon vertex. 

The WTI for the fermion-photon vertex is 



(pi -P2) u T» (p x ,p 2 ) - S- 1 (p0 - S- 1 (pa) 



(IV.l) 



This equation constrains only the longitudinal component of the proper vertex and therefore provides no information 
on the transverse part of (pi,p 2 ). 

Ball and Chiu [3jj have given the form of the most general fermion-photon vertex that satisfies the WTI. It consists 
of a longitudinally-constrained part given by 



1 



1 



Z(px) Z(p 2 



m{pi) _ m (pa) 
Z(px) Z(p 2 ) 



(pi +P2) 



1 



1 



Z(pi) Z(p 2 ) 



(pi +P2Y 
pi -pi 



fa) 



p\ 



Pi 



(IV.2) 



and a transverse part which is described in term of a basis of eight transverse vectors Tj 1 (pi,P2) given in appendix 
lAl The full quark-photon vertex can be written as 



8 

p" (p x ,p 2 ) = r£ c (px,p 2 ) K (pi,p 2 ) t? (pi,p 2 ) 

i—l 



(IV.3) 



where V, (pi,p 2 ) are not constrained scalar functions of p\ and p 2 with the correct C, P, T invariance properties. 

In our scheme the coupling of photons to quarks arises automatically by implementing gauge invariance in our 
lagrangian. The usual way to proceed is to introduce path ordered exponentials in the definition of the currents and 



FIG. 3: Four quark- photon vertex F { T 1 (pi,j>2;P3,P4) where i = S,5,p recalls the current from which the vertex originates 



therefore they become 



Js(x) 



d^yGo (y)ip[x + -y 



ip X 



(IV.4) 



J 5 (a;) 



d 4 yG (y)iP[x + -y) Vie 



f/" +llV dz^A^z) 



iT75 V e 



ip\x - -y 



(IV.5) 



J p (x) 



d A yG p (y)- 



r[x+-y 



1 



V 



iiP\x+-y)p V 



*P[x--y 



ip ( x - -y 



(IV.6) 



where the quark charge is Q — e (f.h + 1/3) /2 with n — (0, 0, 1) , pip (x) = $ip (x) + iQfi (x) ip (x) and ip i x ) p = 
d^ip (x) 7^ — iijj (x)fi (x) Q. In this way Js and J p and J 5 2 become invariant under local gauge transformations. 

The evaluation of the integrals in Eqs. I|IV.4IIV.6|I implies a choice of path. The path dependence is implicit 
in the non locality of the interaction and cannot be avoided. The difference of the contribution between two paths is 
a gauge invariant quantity which is associated with the magnetic flux through any closed surface defined by the two 
paths. However, when the photon momentum vanishes the path dependence disappears. 

The quantization of the photon field in Eas. lUV.4l - HV.6(l leads to the new vertex shown in Fig. [31 whose contribution 
has been fully worked out in Appendix ^ The full quark-photon vertex can be constructed in two steps. The first 
consists in the renormalization of the bare quark-photon vertex by the 4 quarks one photon vertex. Let us call the 
new vertex, shown in Fig. 0] Tq (pi,P2). Using Eas. (|A.7|) and (|A.10|) of Appendix lAl we get 

rft ( Pl , P2 ) = ^ - a [(p! + p 2 f V 0a ip, k) + k» Y ob (P, k)] - 



ap tl±tl + p2 f Ypa (p, k) + ft" W pb (p, k)} 



(IV.7) 



with k = P2 — Pi, P = Pl ~^ P2 and Vo Q , Vot, V pa and Y p b given in Eqs (|A.6|) and (|A.11|I of Appendix 1X1 Both the scalar 
and momentum currents contribute to this vertex. 

The second step consists in the insertion of Tg (pi,P2) in the equation for the quark- photon vertex, represented in 
Fig. \5\ which produces the dressed quark-photon vertex, 



iV 1 ( Pl ,p 2 ) = (pi,p 2 ) + «2.g G 



Pi +P2 



-i2g p 



i>i + h 



G,, 



Pi +P2 



2 

d 4 p 



d A p 

jG p (p) Tr 



Go (p) Tr 



(IV.8) 
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FIG. 5: Diagramatic representation of eq. Hi V.8ft . 

It is easy to show that the solution of equation (|IV.8(1 is just 

r"(p llP2 ) = Y»(p x ,p 2 ) . (IV.9) 
In agreement with our previous discussion, the quark-photon vertex can be rewritten in the form given by equation 

CEHJi, 

^ (p! , P2 ) = T% c (pi , p 2 ) + Vl (pi , P 2 ) Tf (p! , P2 ) + V 2 (pi , p 2 ) I* (p! , p 2 ) , (IV. 1 0) 

where 

Vi (pi,p 2 ) = ^^V 06 (p, fc) , (IV.ll) 

P2 "Pi 

V 2 (pi,p 2 ) = ^4V pfc (p,fe) . (IV. 12) 

P2 Pi 

In summary, we have constructed a dressed quark-photon vertex which has all the desired properties. It is gauge 
invariant in the local sense, without being obtained from the WTI. Two new terms appear from the restoration of the 
local gauge symmetry. The simplest expressions for these terms can be built assuming that a straight line joins the 
two points characterizing the non local currents. For that simple input we obtain 

Vi (pi,j?2) 
V 2 (pi,Pa) 

which can be directly evaluated from the form of the quark propagator. 

V. THE AXIAL CURRENT AND f v . 

Lets us now turn to the coupling of the axial current to quarks. The WTI in this case is 

(P2 [VUP2)T 1 = S- 1 (pajTBT* + 75 r i ,S- 1 (pi) - 2m,^ r 5 (pi^T* . (V.l) 

The last term, proportional to mo, disappears in the chiral limit. 



p\ -pi J -1 



-1 



2 2 



d\ A 



d\ A 



d m (p) 



dp 2 Z (p) 

d 1 

dp 2 Z (p) 



P 2 =(p-#fe) 2 
p 2 =(p— -jfc) 2 



(IV. 13) 
(IV. 14) 
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FIG. 6: Diagramatic representation of the pion decay process corresponding to eq. 



Our first step is to obtain T 5 (pi,p 2 ) T% , which corresponds to the dressing of the vertex i^T 1 . Now the external 
probe is a pseudoscalar-isovector therefore we will have contribution from the pseudoscalar current, J 5 , 



T5 (Pi,p 2 ) t 1 = i-l$T l + i2g Q G Q (p) ij 5 T J 



d 4 p 
(2n) 4 



G (p) (-)Tr 



' s ! P - g ) ^st- 7 iS (p + -jT 5 (p-~,p+ ~ i r 



with pi = p — § and P2 — p + § ■ This equation can be easily solved obtaining 



(V.2) 



T5 (Pi,P 2 ) = *75 



1 + 2 5o G (p) 



Fo (fc 2 ) 



1 - 2 5o n PS (fc 2 



(V.3) 



where F (fc 2 ) is defined by 
5 ij F (fc 2 ) 



^G (p) Tr tst** S (p - ^fc) t Ts^'i s(p+\k 



(V.4) 



We will obtain the dressed vertex Tg (pi,P2) applying to the bosons the same procedure developed for photons 
in the previous section. The explicit calculation is given in Appendix [5] The final result, given in equation (|B.14f) . 
can be rewritten in the following way: 



It(pi,p 2 ) = fg (pi,p 2 ) 



M (Pi,P 2 ) (pi,p 2 ) 75 + V 2 (pi,p 2 ) Tf (pi,p 2 ) 75 



where 



•Al (Pl,P 2 ) = — 2 

P 2 -Pi 



2G (p) - G (pQ - G (p 2 ) 
(Pi -p 2 ) 2 



with Aob (p, fc) given in eq (|A.9b() . and 



It (Pi,P 2 ) = 2 



1 



1 



Z(pi) Z (p 2 ) 
rn(pi) m(p 2 ) 



1 



1 



+ A ob (p, fc) 



(Pi+P 2 ) A 



(V.5) 



(V.6) 



Z(pi) Z(p 2 ) 



1 

" Z(p 2 ) Z(pi)_ 

(p 2 -pl) M .„ r , , (p 2 -Pl)^ 

2 75 - «2m r 5 (pi,p 2 ) 



P 2 ~Pi 



L+ ^ 2 )75 



(P 2 -Pi)' 



(P 2 - Pi)' 



(V.7) 



is the longitudinally constrained part of (pi,p 2 ) . 

Let us now to look for the pion decay constant. In this calculation only the longitudinal current is needed and so 
we do not need to choose a particular path. There are several equivalent ways for obtaining the pion decay constant, 
depending where the quark-quark interaction is included. The simplest one is through the diagram depicted in figure 
El in which all the qq bubbles are included in the Bethe-Salpeter pion amplitude. The axial vertex includes those 
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vertex corrections which can not be included in the bound state amplitude, which in our case corresponds to those 
depicted in Fig0] Therefore, the pion decay constant is defined by, 



d A p 



(-)Tr 



(V.8) 



where (pi,P2) is given in equation (|B.11|) of Appendix iBl and 0^ (p, P) given by equation l|III.3|) . A direct 
calculation gives the result, 



f* = {\ (1 " 2 5onP5 (P 2 )) Pi (P 2 ) + moPo (P 2 )} 



with Pq given in equation 



and 



Pi P 2 = - 



d A p / / P 



Gn 



P-f ) )T.-[,\S»] 



(V.9) 



(V.10) 



We have two clearly distinguishable situations. In the chiral limit /„. is determined by the term with Pi (P 2 ) ■ In 
this limit, Pi (0) is the integral evaluated in equation (|11.12|) . Using in the latter equations 1111. 5fl and H111.6[l and 
substituting, we obtain 



m(0) 
Z(0) 



(V.ll) 



The last form of /„. is the Goldberger-Treiman relation. The crucial point for obtaining this result is the chiral 
symmetric structure of the interaction, which connects the kernel present in the calculation of f„, associated with 
the pseudoscalar-isovector current present in our lagrangian, with the one present in the evaluation of the mass term 
771 (0) , associated to the scalar-isoscalar current. 

If we work with physical pions (mo ^ 0), the pion decay constant is given by 



f m ° Z? ( 2\ 



(V.12) 



In this case the surviving contribution arises from the pion pole present in T§ (pi,p2). 

Using equations (|V.4(1 , 1|II.15(1 and HI.26() we can obtain approximative expressions for Fg 
to the Gell-Mann-Oakes-Renner relation 

fl = -N f { q q }!^ + O(7n , m l) . 



i*A which lead directly 



(V.13) 



We have seen that we have different expressions for f w in the chiral limit and in the physical case. They originate 
from different terms of the axial current. Nevertheless, chiral symmetry guarantees that l|V.ll|) is the limit of equation 
(|V.12(1 when mo goes to zero. 

In table d we give numerical values for for scenarios SI and S2 previously discussed. They are in reasonably 
good agreement with the experimental results. Moreover, this two scenarios are describing the same physics, even if 
they have a very different origin. 

From reference |42j | we can infer that the pion decay constant is determined using equation (|V.8|I with the approx- 
imated expression for the axial vertex, 



r 5 M (Pi,P2) = r 



1 



Z(pi) Z(p 2 ) 



7 M 75 + I 



1 



1 



Z(pi) Z(p 2 ) 



(Pi 



Pi 



— (Pl+h)l5 

-P2 



0(pi-Pa) 



(V.14) 



In this expression O {p\ — pi) implies that (pi — P2)^ ^5 {pi,P2) is determined up to terms of order (pi — p 2 ) — P ■ 
This expression is the minimal generalization of the bare 7^75 vertex in the case where Z (p) ^ 1. It is straightforward 
to proof that the /„■ obtained in this way differs from the exact one by corrections of order m\. Therefore, equation 
(|V.14(1 gives the right result in the chiral limit and a good approximation to the exact value in the physical case. The 
Gell-Mann-Oakes-Renner relation is also well reproduced in this approximation. Nevertheless, equation l|V.14(l is not 
a good expression for the axial vertex, for instance the pion pole is not present in this expression. So we conclude that 
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FIG. 7: Triangle diagram contributing to the pion form factor. 



^5 (PiiPa) is a good approximation of the longitudinal part of Tg (pi,P2) in the vicinity of the pion mass (P 2 ~ to 2 ), 
as it can be see from equation (|B. 15fl . 

Let us now consider the coupling of the axial current to a line of quarks through the vertex, Tg (pi,P2) ■ This 
coupling contains the direct coupling and the pion pole contribution. We observe from equations 1)V.7|1 and l|V.3(l 
that the pion pole contribution goes trough the (j>\,P2) term when the chiral symmetry is explicitly broken. In 
the chiral limit, the pion pole is present in the term with quark masses in Fg (jpi,p2). The full Tg (pi,P2) has some 
dependence on the choice of path. Nevertheless, the overall procedure preserves all the symmetries and, in particular, 
gauge symmetry. 

The longitudinal part of the axial current, Fg (pi,pa) j is path independent. Using equations l|V.12|) and (|V.13(I . 
and assuming that P 2 < m 2 , it can be written in terms of the pion wave function as 



T l 1 
r 5 (PUP*) T7 = o 



1 



Z( P 1 



m (pi 



1 



Z( P 2) 

to (p 2 



a T " 1 



1 



mo 



Z( P i) 



(j>2 - Pi) 



Z(p 2 
<Pn (P) 

P 2 — m 2 



{P2 - Pi) 



T" 



75 T 



(P2~PlY 
(P2 ~PlY 



1 



Z{ P 2) Z( P1 ) 



(pi +P2Y 



-Pi 



T 



(V.15) 



with pi 2 = p ± hP. Equation l|V.15(l manifests the pion pole explicitly, as is predicted by PCAC. 



VI. ELECTROMAGNETIC PION FORM FACTOR. 



We begin by considering the triangle diagram of Fig. [7| In order to hx ideas let us consider the interaction of a 
photon with a 7r + . With the momenta defined as in the figure we have 



!e2WW (fc 2 ) = J 



d 4 P 
(2^) 4 



Tr 



ilbT + (j}^ [p- -k) iS [p- -P 



75 (r + ) f fa [p+^k] % S[p+^P+^k) iQT» ( p + ^P - ^k,p + tp + t k \ lS [ P+ tp-L k 



d 4 p 
(2^ 



(-) Tr 



1 



1 



1 



i75T + 4>„ [p+ -k) iS [p- -P+-k) iQT^ 1 [p - -P - -k,p- -P + -k 



1 



1 



1 



1 



1 - 1 



iS(p--P--k) z 75 (r+y fa[p--k) iS[p+-P 



(VI.1) 



where fa is the pion Bcthc-Salpeter amplitude given in equation (|III.3(I and T M {pi,P2) is the dressed quark-photon 
coupling given in equation IjlV.lOJI . 

There is another contribution to the form factor to be added to the previous one. The J 5 2 (x) term produces a four 



quark photon vertex evaluated in Appendix^] and given by equation 



This four quark photon vertex allows for 
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FIG. 8: Contribution to the form factor coming from the four quark-photon vertex displayed in Fig. |3] 



an additional diagram shown in Fig. |SJ This contribution is 

d 4 p f d 4 p' 



-i(4«7) 



iS(p'- ip'J i 75 (r+) t fa (p') is(p' + \p' 



iS (p+ -P) i 7s r + fa(p)iS (p- -P 



(VI.2) 



06 



Therefore, the full electromagnetic form factor is F (k 2 ) = Ft 2 ^ (k 2 ) + F^ (k 2 ). 

Global gauge invariance guarantees the right normalization for the form factor, F (k 2 = 0) = 1. For a general 
case, the right normalization of the form factor is assured by the Ward identity, equation l|IV.l|) . provided that fa is 

and we take into account the two contributions arising from Eqs. (fvTTTl and (|VF2)l . 

(fill- f>|) . is independent on the total pion momentum, as it is for our 



normalized by equation 

When the kernel used in the pion BSE, equation 
models, the contribution arising from the triangle diagram, equation (|VI.lf) will assure the correct normalization of 
the form factor. Therefore, in our case this property is guaranteed for the full {p\,P2) vertex given by equation 
(|IV.10|I and for the Ball-Chiu vertex, Tg C (pi,p2) , given by equation i|IV.2|) . In our expression for r M (pi,P2) there 
are additional terms besides (pi,P2) which arise from local gauge invariance. They give contributions to F (fc 2 j 
for k 2 ^ without modifying the value at k 2 = 0. 

Due to the separable nature in p and p' and P -independence of the interacting terms in our Lagrangian, one of 
the integrals in equation l|VI.2|l can be done in a trivial manner generating a pion-2 quark-photon vertex. However, 
when performing the remaining integrals we realize that this diagram does not contribute to the form factor. This is 
not a general result, just a consequence of our particular models. In fact this term will contribute in a significant way 
to the parton distribution |46j| . Its contribution is crucial to guarantee isospin symmetry in the parton distributions 
and to restore the momentum sum rule. 

We now proceed to a numerical comparison between the calculations, using the Ball-Chiu vertex, and our full locally 
gauge invariant vertex, equation l|IV.10J) . (|IV.13(I and (|IV.14|I . We show in Fig. ED the form factor for the two vertices 
in scenarios SI and S2, together with the experimental results ^3,E3- We find no important differences for SI when 
comparing the Ball-Chiu prescription to the full vertex. For S2 the correction, which is small for small k 2 , becomes 
important for k 2 ~ — 0.8GeV 2 . The difference in the calculations arises because in S2 Z (p) ^ 1. The correction due 
to Vi (pi,P2) is about 5-7% at this momentum transfer for both scenarios. The correction due to V2 (pi,P2), present 
only in S2, is about 24% for k 2 ~ — 0.8GeV 2 . Since the two corrections go in the same direction the overall result 
changes by about 30%. We have confirmed this conclusion introducing a Z (p) different from 1 in case SI. 

Regarding the experimental results we observe that the scenario SI reproduces well the value of F (fc 2 ) for small 
values of k but underestimates the form factor for k 2 ~ 0.8 GeV 2 by as much as 12%. For scenario S2 we observe that 
the introduction of the full vertex produces a better description of the form factor for k 2 ~ 0.8 GeV 2 but a worst in 
the small k region. In this last scenario the difference between the calculated form factor and the experimental data 
is always less than 5%. 

The behavior for small k can be analyzed in terms of the pion radius. In table [I] we give the mean squared radius 
for the full electromagnetic vertex and, between brackets, that of the Ball-Chiu prescription. We observe that the 
radius is smaller than the experimental result in all cases. We also observe that the full vertex and the Ball-Chiu 
prescription produce the same values since there is no wave function renormalization in the quark propagator, as 
in the case SI. In model S2, with non vanishing wave function renormalization, we observe a difference, of about 
15%, due to the V2 (pi,P2) term in equation l|IV.10|) . confirming our conclusion from the analysis of the form factors. 
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FIG. 9: Comparison of the Pion form factor calculated in the two defined scenarios with the Ball-Chiu ansatz and the full 
vertex of the model. Dot-dashed curve corresponds to the scenario SI with Ball-Chiu ansatz; full curve corresponds to the 
same scenario with the full vertex. Dotted curve corresponds to S2 with the Ball-Chiu ansatz; dashed curve corresponds to S2 
with the full vertex. Experimental data have been taken from |4s| (points) and from |47| (circles). 

Summarizing, the use of the full vertex increases the differences between the calculation and the observation. But 
this is not unexpected since no vector mesons have been included in the models and previous work indicates that this 
contribution can be of the order of 10-20% [HEHi^. 

We analyze the dependence of the pion radius in m (p) and Z (p) in the chiral limit (mo = 0) . In Fig. we show 
the result of rescaling by a generic factor of A the parameters A m and a m appearing in Eas. ijII.24f) and (|II.25[l . If 
A m increases the interaction J 5 2 (x) becomes of shorter range. If we increase a m , go increases. In both cases the pion 
becomes more bound and its radius smaller. 

In Fig. I1UI we also show what happens to the mean square radius when A z and a z are rescaled by a factor of A. 
The system is not very sensitive to changes of A z , while it is quite sensitive for the rescaling of a z when A ~ 2. The 
reason for this strong effect is that for this value Z (0) = 0. 



In the previous section we have analyzed some numerical aspects of the pion form factor. We have put special 
emphasis in the discussion of the terms restoring the local gauge symmetry, Vi (pijPa) an d V2 (pi,P2)- We have seen 
that we have a second contribution given in equation (|VI.2f> . but this contribution vanishes in our particular model. 
Searching for an observable which is sensible to this term, we next discuss the parton distribution. As it is shown 
in |5l| the operator for the parton distribution can be connected to the electromagnetic operator at zero momentum 
transfer. Therefore we are dealing with a property for fc 2 = 0, where the path dependence is absent and our results 
will be valid for any model. 

The first step is to obtain the electromagnetic operator at zero momentum transfer. This operator has a one body 
term which can be obtained directly from Ward identity 



VII. PARTON DISTRIBUTION. 



r*(p,p) 



OS- 1 ( 
dpn 



(P) 



(VII.l) 



and a two body term which corresponds in our particular lagrangian to the equation 1VI.2|) . To be precise, let us 
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FIG. 10: Sensitivity of (r 2 ) in fm 2 for the pion in relation with the propagator parameters in the scenario S2. On the left we 
have lr 2 \ in relation with A • A m for A = 0.5, 1, 1.5 and 2 (dashed curve) and (r 2 ) in relation with A ■ a m for the same values of 
A (full curve). On the right, the dotted curve represents the (r 2 ) in relation with A ■ A z and the full curve gives the (r 2 ) for 
the pion in relation with A ■ a z , for A = 0.5, 1, 1.5 and 2. The full curve gives the msr in fm 2 for the pion in relation with A ■ a z 
for the same values of A. 



proceed to a general discussion. We start from an action of the form 

d^X! d i x 2 d 4 x s d^Xi {xi, %2, £3, x±) (xa) ipp (x 2 ) V> 7 (a; 3 ) tp a (xi) , (VII.2) 



where the greek indices characterize all symmetries, i.e spinor, color and flavor. Let us define the following variables, 

X = (X3 ~ X\) , x' = (X4 - X 2 ) , 

X' — - (x 3 + X\ - x 4 - x 2 ) , X = - (x 3 + X\ +x 4 + x 2 ) ■ (VII.3) 
Translational invariance imposes that Gap-ys (xi, x 2 , x 3 , X4) cannot depend on X. We introduce 

G a p lS {x u x 2 ,x 3 ,x i ) = G ct p l5 {x,x , ,X')= f d \ dAp ' e- tx Pe- lx 'p'e- lX ' p G a07S ( P ,p',P) . (VII.4) 

J (2tt) (27r) (27t) 

Hermiticity and all internal symmetries, such as parity, charge conjugation and time reversal, impose relations between 
the different components of the interaction term G a /3y5 ip,p' , P) ■ 

We are interested in the mesonic bound state described in Fig. [3 The Bethe-Salpeter amplitude for this meson is 
defined by (we identify p 2 = p' + P/2, p A =p' - P/2, p 3 =p + P/2, px=p- P/2) 

r£(p,P) = -2i J ^G a w(p,p',P)(is(p> + ip) T M ip',P)is(p'-\py^^ , (VII.5) 

with the standard normalization condition given in the Appendix IO equation (|C.1|I . 

Regarding the charges of the fields in equation ljVII.2(l . we can assume that Q a = Qs = Qi and Q 7 = Qp = Q 2 . We 
observe that the interaction term in equation i|VII.2(l is invariant under global gauge transformations but not under 
local gauge transformations. One way to make it locally invariant is by incorporating some links. The interaction 
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term becomes, 
S, 



d 4 x d 4 x' d 4 X' d 4 X G aPlS [x, x', X') ip s ( -x' - -X' + X ) e~ lQl I- 



iz^A^z) 



2 2 



2 2 1 



dz^ A^ (z) 



(VII.6) 



We now expand this expression in powers of the photon field. The first term is already included in equation l|VII.2fl . 
The second term is linear in the photon field and is the one of interest. We can evaluate it in the k — > limit without 
having to define a specific path for the integrals present in equation l|VII.6p . From this last result it is easy to obtain 
that the quantity to add in the limit of k — > to each vertex of the type of Fig. [31 is 



Ql ^l,^ (Pl,P3]P2,Pi) -iQt ^2^' (Pl,P3]P2,Pi 



a/3, 7(5 



(4qry) 



a/3, 7(5 



(VII.7) 



with 



r i 4 u 7) (Pl,P3',P2,P4:) 



a/3, 7(5 



dpi dp% 



G 



aft-yS 



Pl + P3 P2+ P4 P2+P3-P1 - Pi 



r 2 4 u 7) {PUP3-,P2,P4) 



a/3, 7(5 



dp% dp% 



G, 



2 ' 2 ' 2 

Pl +P3 P2+ Pi P2+P3-P1- Pi 



(Vlli 



The details of the calculation are given in Appendix [U] where we also discuss the implication of this term to the form 
factor for a BSE with a P— dependent kernel. 

In |5lj . the operator for the parton distribution has been connected with the electromagnetic operator in the 
following way 

d 4 p 



q(x) 



(2k)* 



■Tr 



1 



iS(p--P) T M (p,P)iS [p+ -P 



1 



r f t[p + -p,p + -p) nH s[p + -p) r M (p, p) 



x --.(2p + P) 



(VII.9) 



where T M (p,p) is the dressed photon vertex of the selected parton, equation (|VII.1|) . This expression can be generalized 
from the one body coupling to include the two body coupling in the following way 



'(*) = 

4' 



d 4 P 1 


d 4 p' 


(2ir) 4 J 


(2nf 


v\p,p 


-¥) 



iS[p' --P)t M (p,P)iS[p' + -P 



a/3, 76 



1 



iS[p+-P\T M {p,P)iSlp--P 



1 



(VII. 10) 



with 



[O u (p uP3 ;p 2 , Pi )} af3l5 = - [T, t ( P2 , P3 ) S(x-^(p2 +p 3 )) H) S£ (pi) (2tt) 4 S 4 (p. 

$ (x - I • (pa + P'+Pi)) +5 (x- I • (p 3 + P + Pi))] [ r 2 4 " 7) (P^Ps;P2,P4) 



4 ~Pi 



a/3, 7(5 



(Villi) 



The corresponding diagrams are those shown of Figs. Illlandll2l 

We can split the parton distribution into two one body and one two body contributions 

q (x) = qW (x) + g« ( x ) + q W ( x ) . 

The first contribution, q^ (x) , is the one appearing in equation l|VII.9|) and corresponds to the diagram shown in 
Fig. ^2 The second term is also a one body term, 



(VII. 12) 



fi«(.) = y^(*-n.tt, + P)) 



(2 



Tr 



dT M (p, P) 
dPn 



iS[p+-P) T M ( P ,P)iS [p- -P 



iS[p-^P)f M (p,P)iS[p+^P)in 



dT (p, P) 
dPn 



(VII. 13) 
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FIG. 11: Diagrammatic contributions to the parton distributions of the pion. 





FIG. 12: Diagrammatic contributions to the parton distributions of the pion. 



while the third term is a genuine two body term given by, 



9 (2) [x) 



d A p f d 4 p' 



12„)*J (27T) 



2 M 



i s Ip' - -p r M (p', p)is(p' + -p 



dP u 



+ S[x 



(V + P) 



d 
dPn 



cry 

G a p 7 5 (p',p,p) 



,S( p +lp\ T M (p,P)is(p-±P 



(VII. 14) 



06 



These last two contributions correspond to the diagram of Fig^]. In (x) the bubble integral has been performed 
using the BSE, that is the reason for appearing as a one body term. 

The whole contribution coming from the diagram of Fig lF2l has a non- vanishing contributions to q (x) , nevertheless, 
when we integrate over x the only non- vanishing contribution that can survive is the one associated to the derivative 
of the total momentum, which contributes to the form factor. An analysis of this operator in the scenarios here 
considered is done in |461 |. 

Summarizing, we have three contributions: the standard one, associated with the handbag diagram, and defined 
in equation (|VII.9|> and two new contributions defined in Eqs. (|VII.13() and i|VII.14(l . These contributions are a 
consequence of the non locality of the currents involved in our model. From the point of view of QCD they are also 
handbag diagrams, but in terms of the BS amplitudes they have a different structure. 



VIII. CONCLUSION. 



We have defined a family of phcnomenological chirally invariant non local lagrangians to describe hadron structure. 
They lead to non trivial momentum dependencies in the quark propagator parametrized by momentum dependent 
quark mass terms and wave function renormalization constants. We have shown that the formalism is able to emulate 
any propagator obtained from more fundamental studies. In particular we have studied two scenarios obtained from 
low energy models of QCD ^3[j an( i lattice QCD 0, 0|. As a first goal, our lagrangian description allows a 
careful study of the properties of any observable. In particular we have applied it in detail to the study of the pion 
electromagnetic form factor with special emphasis to the consequences of the local gauge invariance. 

Several authors have previously used non local models, in particular the non local generalization of the Nambu- Jona 
Lasinio model 0, I2H 15^ . l53| and the instanton liquid model [54[ , and have studied electromagnetic properties 
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49, 55]. There are main differences between these approaches and ours. The momentum dependence in our case 
arises from QCD and appears not only in the mass but also in the wave function renormalization. Moreover, we do 
not use a separable approximation of the interaction kernel in each particle momentum. 

In order to test the consistency of the model, we have studied the basic properties of the theory, i. e. effective 
current masses and quark condensate. We have constructed, using our models, the two body quark anti-quark bound 
state equation and solved for the pion obtaining its mass and its Bethe-Salpeter amplitude. 

We have implemented the electromagnetic coupling in this theory. The construction of the dressed quark-photon 
vertex is complicated and the WTI do not solve the problem completely. In particular the transverse propagator has 
been a subject of much debate. In our lagrangian formalism it is natural to implement local gauge invariance through 
links between the points where the quark fields act. We find that two new terms appear in the quark-photon vertex 
which restore local gauge symmetry. We have obtained simple expressions for these two new terms choosing the link 
between the two points characterizing the non local currents as a straight line. They appear as derivatives of the 
mass, m (p), and wave function renormalization, Z (p), of the quarks. We apply our formalism to the pion form factor 
and find out that the local gauge restoring terms could amount to as much as 20-30%. Our analysis is applicable to 
previous work, where the Ball-Chiu prescription for the electromagnetic vertex were used. 

We have applied the same ideas to the axial current, implementing local gauge invariance under SUl (2) transfor- 
mations in the lagrangian. We have obtained the full dressed axial vertex. We have discussed several equivalent ways 
for calculating the pion decay constant. The calculated pion decay constant is in reasonably good agreement with the 
experimental result. Moreover, we can observe that the two scenarios studied in the paper describe the same physics. 
It must be emphasized that this two scenarios have a quite different origin. 

The basic relations from chiral symmetry, the Goldberger-Treiman and the Gell-Mann-Oakes-Renner relations, have 
been recovered. 

The simplest expression for the axial current, given in equation ljV.14fl . is a good approximation for the calculation 
of /jr. Nevertheless it does not include the pion pole contribution. Equation (|V.15|> provides an expression for the 
longitudinal part of the axial current consistent with PCAC. 

Another effect of local gauge invariance is that a new vertex with four quark lines and one photon line appears. 
The contribution of this vertex to the pion form factor vanishes for our models. This is not a general result, but a 
consequence of the separable nature and P— independence of the interacting terms in our lagrangian. For a general 
kernel we have seen that this term is necessary in order to guarantee the normalization of the form factor, F (0) = 1. 

Searching for an observable which is sensitive to this term we have studied the parton distribution. Following the 
ideas of ref [HJ we have obtained the related operator which is path independent and therefore can be used in any 
other model. The new two body term will give a non vanishing contribution to the parton distribution even within 
our models |46j . 

Regarding the numerical values we observe that the two studied scenarios give quite similar results. We only fit mo 
in order to reproduce m^. Their analytic form for the mass, given in equation (|II.23I) and ljll.24f> . are quite different, 
but their numerical value in the region p 2 E < 6 GeV 2 is similar. We regard these expressions as approximations of 
more realistic expressions for the propagator. In this way, their unwanted analytic properties are not considered. 
We conclude that both scenarios provide an overall good agreement with the data. The value of the condensate is 
well reproduced, the value of uiq is in agreement with the quark current mass, the value of the pion decay constant 
is 10% smaller than the experimental value, and both scenarios give similar results. We observe that the fact that 
Z (p) ^ 1 in the second scenario is not important in these observables. The mean square radius discriminates between 
these scenarios. In the SI (S2) we obtain a value which is 10% (20%) smaller than the experimental one. Our 
analysis shows that this difference is not associated to the different choice of the mass expression but to the fact that 
in the S2 Z (p) ^ 1. The electromagnetic radius of the pion will be affected by the coupling of the photon to the 
meson vectors. This contribution has been estimated to be around the 10-20% |38L IHoj. The intermediate region, 
k 2 ~ 0.6 — 0.8 GeV 2 , will be also affected by these vector currents through axial components in the pion Bethe-Salpeter 
amplitude. Therefore we cannot conclude which of the two models is more accurate at present . 

In summary we have set the formalism for the description of non local models of hadron structure and have used 
them to analyze past developments and propose future studies. In particular, it has been already applied for the 
study of the parton distribution of the pion [46|. 
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APPENDIX A: THE QUARK PHOTON VERTEX 



Let us detail here some intermediate steps related with section HVI 

The set of eight basic transverse tensors (p\,P2) we use m equation l|IV.3|) are those defined by Ball and Chiu 
|39| with some minor modifications (k = p 2 — Pi), 

Ti (pi,P2) = Pi (p2-k) - P2 (pi-k) , 

T£ (PUP2) = \p1 (P2-k) - p% (pi-k)} — ^ ^ , 

T£ (pi,p 2 ) = ~i [Pi (P2-k) -p% {px-k)}piPl(T\ v , 

Tg fa,p») = -io^kx , (A.l) 
(pi, pa) = Y {pi - p\) - (pi + p 2 f ft , 

Tj (pi , P2) = \r {h pi 1 - P2 1 } -i(pi+ P2T pIpIvxv , 

Tg (pi,P2) = vy^PxP&xu, + p^i>2 -p^h ■ 

Looking at the four quark photon vertex, we must expand the exponential with the (z) field in the currents of 
Eas ily. 411 V.6|) . Let us start with the scalar current equation (|lV.4fl 

J S (x) = 4 M) (x) + 4 A) (x) + .... - / d^yGo (y) j(x+ ly) $(x- \y\ + (A.2a) 



2 y V 2 



d yG {y)tj} ( x+ -y 



-iQ 



dz^A^ (z)\ ip[x--y 



(A.2b) 



Inserting this expansion in the Jg (x) J$ (x) term of equation (|II.4|> . the first crossed term, Jg (x) J s ' (x) + 

J^ A ^ (x) J^ A ^ (x) , leads to the 4 quarks-photon vertex. An evaluation of this vertex needs to specify the path from 
x — to x + |t/ followed by the evaluation of the integral on 7J 1 in equation l|A.2b|) . This path can be parameterized 
as 



\y^ + cl(X) with AG [-1,1] and eg (-1) = eg (1) = 



(A.3) 



The current J s lA ^ (x) between states of 1 quark of momentum p 1; and 1 photon-1 quark of momentum k and P3, 



respectively gives 



(p 3 \4 1A) (*)| Pifc) = u (p 3 ) (-iQ) u (px) e *(P»-*-Px)«e^ (fc, J .7) 



dH 



u (p 3 ) (-iQ) u (pi) e*^-*-' 1 )*^ (A, 



e ikcs 



P1+P3- Afc 



iC (t) + c^Go (t) 



(A.4) 



t — 



where Q = e (r.n + 1/3) /2 is the charge of the quark, G' (t) = dGo (t) /dt 2 and e M (&,£) is the photon polarization. 
For the sake of simplicity we take cs M = which corresponds to linking the two points of the non local current by a 
straight line. We can write the matrix element in the following way 



(p 3 \4 1A) (x)\ Pl k) = -ti (pa) Q u (pi) eKw-fc-pOx^ M 



(Pi +P3) M V a ( Pl P3 , k I - /,■„ 



Pi + P3 



■ A' 



(A.5) 
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where 



V a (P, fc) 



(A.6a) 



V„b(p,*) = / ^A-G (p-A- 



(A.6b) 



Using the current given in equation 1A.5JI we obtain the vertex associated to the interacting term 
.go (4° A)t («) 4 1A) (*) + 4 1A)t («) 4° A) (*)) of the Lagrangian, 



= -«2.g Go ( P2 \ ) SpsQya 
Q/3,7<5 V 2 



(Pi + P3) u V 0cl ( Pl - P3 , fc 



kuVob 



Pl +P3 



+ (Pia,P37 P2fl 1 Pi&) , 



(A.7) 



with fc = pa + P4 — pi — P2- 



We apply the same ideas to the pseudoscalar current i|IV.5(l obtaining the vertex associated to the interacting term 
.90 (4° A)t (&) 4 1A) (*) + 4 1A)t ("0 4° A) (a:)) of the Lagrangian, 



r 5 4 u 7) (Pl,P35P2,P4) 



a/3, 7(5 



-i2 3o G ( P2 n ^ 4 ) (i75T) 5i g 



*75 2 W^l 



7CK 



(Pl+flB) M Vo.(^,fc)-fc„V 06 f^±^.* 



JOC 



Pl ~\~ P3 

(Pl +P3) U A 0a ( — ^ — ) - k^Aob 



2 

Pl + P3 



+ (pi«,P37 ^P2(3,pa5) 



with 



(P, *) = jf I G[) (p - A|) - £ dA i Go (p - A^ 



(A.8) 



(A.9a) 



A06 (P, fc) = jT 1 dA A Go ^ _ A ^ - 1° dX ~G' Q (p- X^ 
and Voa and Vot given in equations l|A.6al) and l|A.6b|) . 



(A.9b) 



The momentum current equation (|iV.6(l produces the interacting term g p (^Jp nA ^ (x) Jp 1A ^ (x) + Jp 1A ^ (x) Jp ^ (x)^j 
The associated vertex is 

'P2 + Pi \ fj>2+l>i 

a/3, 76 



whith 



r l 4 u 7) (Pi,P3;P2,P4) 



(Of) 



70; 



9*n 2 



Pi + P3 + k 



Pl + P3 - k 



(Pl +P3) V 



(Pl + P3) u Vpa ( Pl t P3 1 fc ) ^ 



" * 2 

Pl + P3 



{Pl(±,P3l ^P2f3,P48) 



V pa (P, k) 



J^dX^G' p (p-X 
Y pb (p,k) = J ^dX±G' p (p-X 



(A.10) 



(A.lla) 



(A.llb) 
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APPENDIX B: THE QUARK W" VERTEX. 

Let us decompose the quark field in its left a right parts, ip L ' R (x) = i (1 =F 75) ip (x) . Under a local SUl (2) gauge 



transformations we have 



iP L (x) — > e ifc ?- 5(l V L (x) ip R (x) — ► (a:) 
tW' 1 (1) = tW^ 1 (x) - <9^a (x) .r - g w (a (x) x (x) ) .f 



(B.l) 
(B.2) 



Our lagrangian model defined in (|II.4fl is invariant under infinitesimal global gauge SUl (2) transformation. Now we 
are interested in making this lagrangian invariant under local transformations. We proceed in a similar way as in the 
electromagnetic case. The main difficulty is on the fact that our currents are built with fields in different points, thus 
the value of a(x) will be different for each point. To avoid this difficulty let us introduce the following product of 
fields 



iSML J* dz* TW^(z) 



r (») 



(B.3) 



which transform covariantly at least for infinitesimal transformations, 



Of. r dz" ?w»(z) , l 



V (y) 



e i2f-T.a(x) 



iM ¥ .J*dz<* TW»(z)+0(gl) ,L 



V (y) 



(B.4) 



Then, we can define the currents 

J s (x) = f d 4 yG (y) 



ip n [x+-y) V[e 



-i^f^idz^W^zy, 



i> ( x - -y 



(B.5) 



,/ 5 (x) = / d 4 yG (y) 



-lip n \x+-y\ rPle J" 1 ip \x--y 



i ip v [ x+-y ) V\i 



dz^W^(z)T 



Tip 1 



x - t;V 



(B.6) 



j p ( x ) = / d 4 yG p (y) ■ 



^ \x+-y) V\e 



-i^lXdz^(z)r\ f 1 

2 " I %P tp I x - -y 



# [x + -y\ p V e J — *» 



( a; + J y V fl U - ^ 



$ ( a; - -y 



(B.7) 



where the covariant derivative is D^ip (x) = <9 M + ig w W^ (z) 5 ip L (x) . With these definitions, the combination 

Jj (x) + J 5 2 (x) is invariant under local infinitesimal gauge transformations. The momentum current, Jp (x) , is 
self-invariant. 

Once we have the invariant lagrangian we proceed for obtaining the 4quarks-M^ /i vertex. The procedure is exactly 
the same as in the previous appendix, expanding the currents in the number of W mesons, and we write directly our 
results 

The vertex associated to the interacting term g ^ J^ W ^ (x) J$ 1W ^ (x) + Jg 1W ^ (x) J^ W ^ (x)\ of the Lagrangian, 
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in which a boson of type W 1 is involved, is 

" r ,(4gW) / \ 
- ^Vob 



a/3, 7(5 

Pl + P3 , 
2 ' 

Pi +P3 



" 2 ^ffO-^-CrO 



P2 +P4 



(Pi +P3)„ V 0a 



(Pi +P3) U A 0a 



Pi +P3 



Pi +P3 



,k 



(B. 



The vertex associated to the interacting term go ( J^ 01 ^ (a;) Js 11 ^' (x) + J^ 1W ^ (x) J^ W ^ (x) ) of the Lagrangian is 



r 5, M (Pi-P3;P2,P4) 
Pi +P3 



2 

Pi + P3 



- % 2ga—G 



P2 +P4 



(H5T 3 ) 



US 



(Pi +P3)„ V 0a 



Pi +P3 



70; 



(Pi +P3) M A 0a 



Pi +P3 



(B.9) 



The vertex associated to the interacting term <? p ^ Jp 0W ^ (x) Jp 1W ^ (x) + Jp 1W ^ (x) Jp 01 ^' (xfj of the Lagrangian is 



ltf°(^3;P2,p 4 ) 



q/3,7<5 

i 



G,, 



P1+P3 + k 



G r , 



rVo(l-7s) 



7a 



Pi + P3 - 



(Pi +P3)^ 



(Pi +P3) il V pa ( Pl + P3 , k 



(Pi",P37 ^P2/?,P4<5) , 



(B.10) 



In all these expressions, the momentum transferred is defined as fc = P3 + P4 — Pi — P2- The functions Vo a , Vo&, A 0a , 
Aob, V pa and V p b are defined in equations (|A.6a|) . l|A.6b|l . l|A.9a|l . I|A 9bfl . I|A. llaf) and (|A.llb() . 

The H^ 1 boson couples to quarks through the dressed vertex (r M (pi,P2) — ^5 (pi,P2)) t 1 . As in the quark photon 
vertex discussed in section llVl the full quark W vertex is constructed in a two steps process. The first one is shown 
in fig. 0]and it consists in the renormalization of the bare quark- W vertex by the 4 quarks one W vertex. The change 



in the vector part of the vertex is given in equation i|IV.7|) . The axial part is 

r 5,o (Pi.Pa) = 7 M 75 + "o [(Pi +P2T A 0a (p, k) + k» A ob (p, k)} 75 + 

'Pl +P2\ (Pi -P2) M 



2 Go 



(pi -P2) 



r75 [«o + goFi (k)] 



ap-[G p ( Pl ) + Gp (pa)] 7^5- 



a P — n — 75 



(pi+p 2 f Vp Q 



Pl +P2 



k » ¥pb[ P±±Pl, k 



with F\ (k) given by equation ijV.lOf) . It is interesting to note that the longitudinal part of T§ (pi,P2), 



(P2-Pi)u r 5,o(Pl>P2) = S 1 (P2)l5 +I5S 1 (pi) + [ 2m ~ 2g G 



Pi +P2 



^i(fc 2 ))75 



(B.ll) 



(B.12) 



does not satisfy the WTI (fVTjl . 

The second step is represented in Fig. [5j The associated equation for the vector part of the vertex is given in 
equation I|1V.8|I and the final form of the solution is given in equation i|lV.10|) . The axial part of the vertex is 
governed by the equation 



T 1 T 1 

(P1.P2) y =^^0(^,^2)—- 



i 2 g G 



Pl + P2 



? 75 T J 



d A p 



G (p) (-)Tr 



i S 



(B.13) 
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with k = P2 — Pi ■ This equation can be solved obtaining 

r£ (pi»P2) = r£ (pi,p 2 ) + 2.90G0 f^liJ? 2 . 

This expression can be rewritten in the form given by equation (|V.5 



F (k 2 ) 



^(fc 2 )+2m \,, 2 
1 - 2.g Hps {k A 



k^ 

F 75 



(B.14) 



In equation l|V.14f) we give an approximated expression for the longitudinal part of Tg (pi,P2) • From equations 
(|V.14fl and (|B.12j) we have 



(pi -k)^ (r£ (pi,p 2 ) - r' 5 f (pi,p 2 )) = 

[ao G (pt) + a Q G (pa) + 2 ffo G (p) F x (P 2 )] 7s 



(B.15) 



with px,2 — p±P/2. Inserting equation (|B.15|I in (|V.8|) we can evaluate the numerical error produced in the calculation 
of f n through the use of r g ^ (pi, P2) • Expanding in powers of P 2 it is straightforward to obtain that A/ ff = O (m 2 ) . 

There are alternatives ways for calculating the pion decay constant in which we cannot use the approximated 
expression (|V.14fl . For instance, we can consider an interacting pair qq which a some point couples to a W boson. 
We can describe the interaction between the qq by the scattering amplitude or using the dressed vertex. This two 
descriptions must be equivalents and in the proximity of the pion pole we have 



(27r)" ^ 2 J 3 V 2 2 7 2 V 2 

Where, in our model, the interaction is 

G aPl s {p',p, P) = go Go (p') G (p) 8 1OI S S + (ifj 5 )^ a (ifj 5 ) S/3 + g v G p (p 1 ) G p (p) (j>') ia (j>) sp 
and the pion amplitude is 

K J (p,P)= i g« qq G a (p) (*7 5 r 4 ) Q7 



06 



(B.16) 



(B.17) 



(B.18) 



In the right hand side of equation (|B.16|) we must consider only the pion pole contribution. It is easy to reproduced 



the result obtained in equation (|V.9(1 for f„ , and is also obvious that we cannot use of the expression (jV.14(l , because 
we lost the pion pole in the right hand side of equation l|B.16|) . 

APPENDIX C: FOUR QUARKS-PHOTON VERTEX IN THE GENERAL CASE AND THE VALUE OF 

THE FORM FACTOR AT k 2 = 0. 



The standard normalization condition for the Bethe-Salpeter amplitude is 
2iP f * = 



d 4 p 
(2^ 



Tr 



t M {p,P)i 



-2i 



dS ( P + |P) 



where 



dP u 



r'VPi l *{ P -±p\+T M (p,p)iSlp + ^p) r- w (p.P) / 



dS( P -\P) 



4^/ r 



<Pp <Pp 
(27T) 4 (2^) 4 



<S'(// + V) T M (p',P) is(p' - l -P 



f M (p,p) = 7o [r M (p,p)] t 7o 



0:7 



06 



dGgpjs (p,p',P) 

OPa 



dP u 



(C.l) 



(C.2) 



This normalization condition is equivalent to eauation l|III.6|) 
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A minimal test of consistency of our calculation is that the form factor at k 2 = must be 1 for an amplitude 
normalized with equation (|C.lf) . Usually the form factor is calculated in the impulse approximation, which includes 
only the triangle diagram shown in Fig. The use of the Ward identity l|VII.l|l in order to define the electromagnetic 
vertex in this diagram provides a contribution which coincides with the first integral on the right hand side of equation 
(|C.1() . From that we can conclude that the use of the BSE for the pion with a P-independent Bethe-Salpeter kernel 
together with the triangle diagram provides a consistent approximation scheme 36]. For a P-dependent kernel this 
consistency is lost even at k 2 — due to the presence of the second integral on the right hand side of equation (|C.1(I . 

Let us proof that F (0) = 1 from equation l|VII.6|) . We need to expand equation l|VH.6[] in powers of the photon 
field. We retain the second term, which is linear in the photon field. We can evaluate it in the limit k 2 — > without 
defining a specific path for the integrals present in eauation (|V11.6ll . obtaining 



Qi 



X 



dz fi A^(z) + Q 2 



dz^Af, (z) + Q 2 
-+x Jx 



-+X r X 

dz^A^ (z) + Qt 



k 2 ^Q 



6p 



dz^A^ (z) 



-x 



(C.3) 



From this last result it is easy to see that the quantity to add to each vertex of the type of Fig. [21 in the limit of 
k 2 — > 0, is given by Eqs. HVII.7|) and IjVII.Sjl. 

The electromagnetic form factor in the limit k 2 — > 0, including the contributions from Figs and [SJ is 



ie2P^P(0) 



-1Q2 



d 4 p 
d 4 p 



Tr 



Tr 



T M (p,P)iS \p--P) T M (p,P) iS (p+ -P [p+ -P,p+ -P) iS [p+ -P 



T M (p,P) iS [p- -P r** [p - -P,p - -P iS[p--P) T M (p,P)iS[p+-P 



d 4 p f d A p' 



(2nTJ (2tt) 



% S [p 1 - -P )T M (p', p)iS[p' + -P 



1 



ay 



-iQi 



^(497) ( „/ 



p' - -p,p' + -p ;p +-p : p--p 



"" (p'-lp,p' + lp ]P + lp, p -lp 



2,fi 



a/3, 75 



iS(p+~p)r M (p,p)iS( P --p 



(C.4) 



In order to simplify the first two lines of this equation (the one body part) we make use of the WTI, equation \VR.l) . 
For the two body part of the equation we use 



iQiTff [ p >-lp, p > + ~P ]p+ ~P,p-lp 



1 



-1Q2T 



a/3, 7(5 



-iQi 
-iQ2 



d 


d 


-2 4 


dp'v 


^ dpi* 


dP^ 


d 


d 


„ d 
+ 2 


dp'v 


h dp^ 


dP^ 



Equation 1VII.5|) allows to do some of these integrals. Charge conjugation symmetry leads to 



j£p_ 



Tr 



d r M (p, p) 



dp» 



iS[p+-P) T M (p,P)iS[p--P 



d 4 P 



Tr 



iS[p--P)t M (p,P)tS[p+^-P 



l\dT M (p,P) 



dpv 



. 



(C.7) 



With these inputs, the normalization condition of the Bethe-Salpeter amplitude, equation (|C.1|) . implies the natural 
normalization for the form factor, F (0) = 1. We observe that only the contribution associated with the derivative of 
the total momentum in Eas. ljC.5fl and (|C.6[l gives a non vanishing result. 

Our results show that consistency between the Bethe-Salpeter normalization condition equation i|C.l|) and the value 
of the meson form factor at k 2 = is also attainable for a P-dependent kernel, if we add the contributions coming 
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from the diagrams of Figs. [7] and [SJ This result is consistent with field theory. We have simply added all the diagrams 
with a one photon coupling, no matter where the photon couples in our system, and in this way we have obtained the 
gauge invariant contribution to the form factor. Fig. [HI confirms that the use of the Ward-Takahashi identities for the 
components of a system is not sufficient to assure that the gauge symmetry is satisfied for the composite system. 
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